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Abstract 

> 

f«\j ' An effect generated by the nonexponential behavior of the survival 

1^ ■ amplitude of an unstable state in the long time region is considered. 

We find that the instantaneous energy of the unstable state for a large 
^^ _ class of models of unstable states tends to the minimal energy of the 

f~^ ■ system Smin as t — )• oo which is much smaller than the energy of this 

^D , state for t of the order of the lifetime of the considered state. Analyzing 

the transition time region between exponential and non-exponential 
form of the survival amplitude we find that the instantaneous energy of 
^ , the considered unstable state can take large values, much larger than 

H I the energy of this state for t from the exponential time region. Taking 

into account results obtained for a model considered, it is hypothesized 
that this purely quantum mechanical effect may be responsible for the 
properties of broad resonances such as a meson as well as having 
astrophysical and cosmological consequences. 
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1 Introduction 

Searching for the properties of unstable states |0) G "H (where "H is the 
Hilbert space of states of the considered system) one analyzes their decay 
law. The decay law, V^{t) of an unstable state \(j)) decaying in vacuum is 
defined as follows 

V4t) = \ait)\\ (1) 

where a{t) is the probability amplitude of finding the system at the time t 
in the initial state \(j)) prepared at time to = 0, 

ait) = (0|0(t)). (2) 

and \(p{t)) is the solution of the Schrodinger equation for the initial condition 

10(0)) = \cpy. 

th^m)) = Hm)), (3) 

where H denotes the total selfadjoint Hamiltonian for the system considered. 
From basic principles of quantum theory it is known that the amplitude 
a(t), and thus the decay law '^^(t) of the unstable state |0), are completely 
determined by the density of the energy distribution u}{S) for the system in 
this state [D 6] 

a{t)= [ u{£)e~'h^td£. (4) 



JSpec.{H) 

where uj{£) > and a(0) = 1. 

In |3] assuming that the spectrum of H must be bounded from below, 
{Spec.{H) > — oo), and using the Paley- Wiener Theorem |1] it was proved 
that in the case of unstable states there must be \a{t)\ > A exp [— fet''] for 
|t| -)■ oo (where A> 0, b > and < g < 1). 

The problem of how to detect possible deviations from the exponential 
form of V(i){t) at the long time region has been attracting the attention of 
physicists since the first theoretical predictions of such an effect. Many tests 
of the decay law performed some time ago did not indicate any deviations 
from the exponential form of 'P,/,(t) at the long time region. Nevertheless, 
conditions leading to the nonexponetial behavior of the amplitude a{t) at long 
times were studied theoretically. Conclusions following from these studies 
were applied successfully in an experiment described in the Rothe paper [5], 
where the experimental evidence of deviations from the the exponential decay 
law at long times was reported. This result gives rise to another problem 



which now becomes important: If (and how) deviations from the exponential 
decay law at long times affect the energy of the unstable state and its decay 
rate at this time region. 

Note that in fact the amplitude a{t) contains information about the decay 
law V^(t) of the state |0), that is about the decay rate 7° of this state, as 
well as the energy S^ of the system in this state. This information can be 
extracted from a(t). Indeed, if |0) is an unstable (a quasi-stationary) state 
then, there is 

h 2^^-'^ dt aoity ^^^ 



for t r^ T^, where 



ao(t) = exp [--^(£:J - ^70)t] ^ a{t), (6) 



for t ~ r^, r^ = 4t and 7^ is the decay rate of \(j)). 

The standard interpretation and understanding of the quantum theory 
and the related construction of our measuring devices are such that detect- 
ing the energy £'^ and decay rate 7° one is sure that the amplitude a{t) has 
the form ([6]) and thus that the relation ([5]) occurs. Taking the above into ac- 
count one can define the "effective Hamiltonian" , h^, for the one-dimensional 
subspace of states l-i\\ spanned by the normalized vector |0) as follows (see, 

eg. i), 

, def .^ da{t) 1 

In general, h^ can depend on time t, /i^ = h^{t). 
It is easy to show that equivalently |7] 

"*'" = 7#Mr- ''' 

One meets effective Hamiltonians of this type when one starts with the time- 
dependent Schrodinger equation ^ for the total state space % and looks 
for the rigorous evolution equation for the distinguished subspace of states 
7^11 d % (see ^ and references one finds therein). In the case of one- 
dimensional "Hi I this rigorous Schrodinger-like evolution equation has the 
following form for the initial condition a(0) = 1, |6], 

'^^^ = Kit) ait). (9) 



Relations ([7]) and (jHl) establish a direct connection between the amplitude 
a{t) for the state |0) and the exact effective Hamiltonian h^{t) governing the 
time evolution in the one-dimensional subspace T-L\\ 3 \(fi)- Thus, the use of 
the evolution equation ([9]) or the relation ([7]) is one of the most effective tools 
for the accurate analysis of the early- as well as the long-time properties of 
the energy and decay rate of a given quasi-stationary state |0(t)). 

So let us assume that we know the amplitude a{t). Then starting with 
this a{t) and using the expression ([7]) one can calculate the effective Hamil- 
tonian h^{t) in a general case for every t. Thus, one finds the following 
expressions for the energy and the decay rate of the system in the state |0) 
under considerations, to be more precise for the instantaneous energy and 
the instantaneous decay rate, (for details see: [7]), 

£^ = S^{t) = ^{h^{t), (10) 

70 = 70W = -2S(/i0(t), (11) 

where 3ft {z) and Q= [z) denote the real and imaginary parts of z respectively. 
Using dZD and ([TU]), (HI]) one can find that 

(12) 
(13) 
(14) 
(15) 

The aim of this talk is to discuss the long time behaviour of S^{t) using 
a(t) calculated for the given density uj{S). We show that ^^^(t) — )■ Smin > — oo 
as t — 7- oo for the model considered and that a wide class of models has similar 
long time properties: £(f,(t)\^_^^ ^ £^. It seems that, in contrast to the stan- 
dard Khalfin effect [3j, in the case of the quasi-stationary states belonging 
to the same class as excited atomic levels, these long time properties of the 
instantaneous energy ^^(t) have a chance to be detected, eg., by analyzing 
the properties of the high energy cosmic rays or the spectra of very distant 
astrophysical objects. 



^0) 


= (0i^0), 


S^{t ~ T^) 


^ £l ^ ^M 


7^(0) 


= 0, 


70(^ ~ T-^) 


^ ll- 



2 The model 

Let us assume that Spec.{H) = [Smin, oo), (where, Smin > — oo), and let us 
choose oj{S) as follows (compare [8]) 

N l" 

u{£) = UBw{£,£min) = 7^ &{£ — £min) o ) (16) 

where A^ is a normalization constant and 

1 for £ > 0, 



^^^^ ~ 'I for £ < 0. 
For such ojbw{.£) using (jl]) one has 

a{t) = ^ r ^° n e-l^t dS, (17) 

where 

T7 = — / —d£. (18) 

N 27rJ^_^^_^0)2 + (^)2 

Formula fITT]) leads to the result (see also |8]) 



— X 

27r 



a(t) = A^e ft*^^<^ S"")^/! 

-i?i(-;^(^?-^™n-^75)t)]}, (19) 

where Ei{x) denotes the integral-exponential function [8l|9]. 
In general one has 

a{t) = aexp{t) + anon{t), (20) 

where 

aexp{t) = A^e~t(^0 - ^^)^, anonit) = a{t) - a.xpit). 



Making use of the asymptotic expansion of Ei (x) [H] 

^i(-)li.i^oo - ^(1 -; + !-•••)' (21) 

where | a,igz\ < fvr, one finds 

N -is ■ t h ■^ 7^ h 



'^'^ ^ \ '''(f) '^mm I '' 



where h^ = £^ - 1 7^, and 



^.(01,-.. = «*<'^ ' 



St a(t) 



t— >oo 



h _ t:-^ t^rjiin / h\'^ 



f' I fA '-'min \ ^ L ' 

for the considered case (IT6|) of w^vkI^) (for details see [3). From ( 123!) it 
follows that 

fO^f. /?^N 2 



E ■ - 1 






i— >oo 



where S^{t) = £^{t)\t^oo, and 

h 



^ {hMt^J ^ - - ^ 0. (25) 



t—^oo 



The property ( 12^ means that 



"^iuti^j = ^rw <^^- (26) 
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For different states \(f)) = \j), {j = 1, 2, 3, . . .) one has 

^ih{t)\t^J = ^ih,itl^J, (27) 

whereas in general 7° 7^ 73. 

Note that from (l22l) one obtains 

7° 

I V /IC^-GO I 

+— sin [(£0 - ^_) t] e-2 - ^ J^ - 

^^^^ ^ + ••• • (28) 



I (/) ^min I t- 

Relations fl22|) — ( 127|) become important for times t > tas, where tas denotes 
the time t at which contributions to |ct(^)|t_^ool from the first exponential 
component in (128!) and from the third component proportional to ^ are 
comparable, that is (see (l20l)). 

|ae..p(t)|2 ~ |a„o„(t)|2 (29) 

for t — 7- 00. So tas can be be found by considering the following relation 

Assuming that the right hand side is equal to the left hand side in the above 
relation one gets a transcendental equation. Exact solutions of such an equa- 
tion can be expressed by means of the Lambert W function (TU]. An asymp- 
totic solution of the equation obtained from the relation (1301) is relatively easy 
to find [11]. The very approximate asymptotic solution, t^s, of this equation 

F F 

for (-5-) > 10 (in general for (-5-) — )■ 00) has the form 

14, T^ 



70 



0+ F^ - F 



h 



8,28 + 4 In 



vO 



+ 2 1n[8,28 + 4 1n( ^ „ '"^" )] + ... . (31) 

7^ 



3 Some generalizations 

To complete the analysis performed in the previous Section let us consider a 
more general case of uj{S) and a(t). For a start, let us consider a relatively 

simple case when limf_+^^.,^+ ^{^) = wq > and uj{S)\£^£ . = 0. Let 
derivatives u^''\£), {k = 0,1,2, ... ,n), be continuous in [Smim oo), (that is 
let for £ > Smin all uj^''\S) be continuous and all the limits lim£-_^£^^^+ u'^^\S) 
exist) and let all these u^^\£) be absolutely integrable functions then (see 

n-l 



a{t) 



ih —if . f 

t 
t— >-oo fc=0 



E(-i) 



Ah, 



k (k) 



CO, 



) 



(32) 



where 4''^ = limf^£„,„+ uj^''\S). 

Let us now consider a more complicated form of the density oj{£). Namely 
let uj{S) be of the form 



u{£) = {£ - £min)^ ri{£) e Li(-oo,oo) 



(33) 



where < A < 1 and it is assumed that r]{£) > 0, ri'^^\£), (A; = 1, 2, ... , n), 
exist and they are continuous in [£min-, oo), and limits lim£-_^£-^^^^+ ri'^^\£) ex- 
ist, lim£:^oo {£ — £minY 'n'^^\£) = for all above mentioned k and 






0, then 



a{t) 



ih , _ip . + 
t 

. 111. 9 , . 

+ (- j) «n-2(t) 

+ (- j) an-3(t) + . . . 



a„(t) + (- — )a„_i(t) 



where (compare [T3l fl^) 

n— fc— 1 



i=0 






T{z) is the Gamma Function and r]^ = lim£;_>£:^.^+ ri^^\£), r]^^\£) 
and j = 0, 1, . . . , n. 



(34) 



(35) 
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The asymptotic form of hu{t) for t — )■ oo for the a{t) given by the relation 
looks as follows 

-^(^^ + .... (36) 

In the more general case of a{t) (see, e.g. (^^ ) after some algebra the 
asymptotic approximation of a{t) can be written as follows 



a(t) ~ e-^^^™"E^, (37) 

t— i>oo k=0 



where ^ > and Ck are complex numbers. 
From the relation fl371) one concludes that 



^^ . e-l^^'^ ■ 



{-^^™"-E(^+^)^}- (3^ 



at 

i— >oo fc=0 

Now let us take into account the relation (|9]). From this relation and 
relations f l37|) . fl38|) it follows that 

Kit) ~ f™„ + ^ + -f + -f + . . . , 39 

t— >oo 

where di,d2,ds, ■ ■ ■ are complex numbers with negative or positive real and 
imaginary parts. This means that in the case of the asymptotic approxima- 
tion to a{t) of the form ( !37l) the following property holds, 

lim h^{t) = S.min < £l (40) 

t— >oo ^ 

It seems to be important that results ( l39l) and ( HOl) coincide with the 
results ( 123|) — ( 1271) obtained for the density a;BVK(^) given by the formula 
(fT6|) . This means that general conclusion obtained for the other u]{£) defining 
unstable states should be similar to those following from (123|1 — ([2] 



4 Numerical calculations 

Long time properties of the survival probability |a(t)p and instantaneous 
energy S^{t) are relatively easy to find analytically for times t ^ tas even 
in the general case as it was shown in previous Section and [12]. It is much 
more difficult to analyze these properties analytically in the transition time 
region where t ^ tas- It can be done numerically for given models (see [15]). 
The model considered in Sec. 2 and defined by the density ubw{S), (IT6|) . 
allows one to find numerically the decay curves and the instantaneous energy 
£(^(t) as a function of time t. The results presented in this Section have been 
obtained assuming for simplicity that the minimal energy Smm appearing in 
the formula fITB]) is equal to zero, Smin = 0. So, all numerical calculations 
were performed for the density UBwi^) given by the following formula 

N Is 

^BW{^) = UJBw{^,£min = 0) = 7— 0{S) -, , (41) 

for some chosen -^. Performing calculations particular attention was paid to 

the form of the probability |a(t) p, i. e. of the decay curve, and of the instan- 
taneous energy e<)i(t) for times t belonging to the most interesting transition 
time-region between exponential and nonexponential parts of |a(t)p, where 
the following relation corresponding with ( 129|) and ( 130|1 takes place, 

\aexp{t)\^ ~ |a„on(t)P, (42) 

where aexpif) ^ o.non{t) are defined by fl20|) . Results are presented graphically 
below in Figs (P) — (CD- 



5 Final remarks. 

Decay curves of a type Fig. ([1]), Fig. ([2]), Fig. (jl]) — Fig. ([6]) one meets for 
a very large class of models defined by energy densities to{£) of the following 
type (see [21 [H]), 
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Figure 1: Survival probability V,f,{t) 
The case 5 = 10. 



|a(t)p in the transition time region. 



where A > 0, f{S) is a form-factor — it is a smooth function going to zero 
as £^ — )■ cxD and it has no threshold and no pole. The asymptotical large time 
behavior of a(t) is due to the term {£—£minY ^"^^ the choice of A (see Sec. 3). 
The density ui{£) defined by the relation f H3|) fulfills all physical requirements 
and it leads to the decay curves having a very similar form at transition times 
region to the decay curves presented above. The characteristic feature of all 
these decay curves is the presence of sharp and frequent oscillations at the 
transition times region (see Figs ([1]), ([2]), (|1]), ([5]), (|6])) (see also, eg. p^ITT]). 
This means that derivatives of the amplitude a{t) may reach extremely large 
values for some times from the transition time region and the modulus of 
these derivatives is much larger than the modulus of a(t), which is very small 
for these times. This explains why in this time region the real and imaginary 
parts of h^{t) = £(f,it) — |7</,(t), which can be expressed by the relation ([7]), 
ie. by a large derivative of a{t) divided by a very small a(t), reach values 
much larger than the energy £9 of the the unstable state measured at times 
for which the decay curve has the exponential form. For the model considered 

we found that, eg. for -n- = 10 and Sr^ < t < GOr,^ the maximal value of 

the instantaneous energy equals (^^(t) = 89, 2209 f^^ and £(f,{t) reaches this 
value for t = tmx, lo = 53, 94 r,^ and then the survival probability V^f^it) is of 
order P^(tm^_io) ~ 10~^. 



11 




Figure 2: Survival probability V(f,{t) = \a{t)\ in the transition time region. 
The case -g- = 10. 



SAt) 




moT(t> 



Figure 3: Instantaneous energy £tf,{t) in the transition time region. The case 
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Figure 4: Survival probability Vt^^t) 
^ = 100. 



|a(f)p in the transition time region. 



The case ^ 




2.x 10 



26 T<^ 



Figure 5: Survival probability V(f,{t) 
The case ^ = 100. 



|a(t)p in the transition time region. 
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Figure 6: Survival probability V^{t) = |a(t)p in the transition time region. 



The case % = 100. 
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Figure 7: Instantaneous energy S^{t) in the transition time region. The case 



4 = 100. 
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The question is whether and where this effect can manifest itself. There 
are two possibihties to observe the above long time properties of unstable 
states: The first one is that one should analyze properties of unstable states 
having not too long values of tas- The second one is finding a possibility to 
observe a suitably large number of events, i.e. unstable particles, created by 
the same source. 

The problem with understanding the properties of broad resonances in 
the scalar sector {a meson problem [19l|) discussed in [201 EI], where the hy- 
pothesis was formulated that this problem could be connected with properties 
of the decay amplitude in the transition time region, seems to be possible 
manifestations of this effect and this problem refers to the first possibility 
mentioned above. There is the problem with determining the mass of broad 
resonances. The measured range of possible mass of a meson is very wide, 
400 - 1200 MeV. So one can not exclude the possibility that the masses of 
some a mesons are measured for times of order their lifetime and some of 
them for times where their instantaneous energy S„{t) is much larger. This 
is exactly the case presented in Fig. (^ and Fig. ([7]). For broad mesons the 
ratio -§■ is relatively small and thus the time tas when the above discussed 
effect occurs appears to be not too long. 

Astrophysical and cosmological processes in which extremely huge num- 
bers of unstable particles are created seem to be another possibility for the 
above discussed effect to become manifest. The probability ^^(t) = |a(t)P 
that an unstable particle, say 0, survives up to time t ~ tas is extremely 
small. Let V^{t) be 

nw !,.,„, - lo-^ (44) 

where k ^ 1, then there is a chance to observe some of particles <p survived 
at t ~ tas only if there is a source creating these particles in J\f^ number such 
that 

nWI*^*..-^^ » 1- (45) 

So if a source exists that creates a flux containing 

Af^ ~ 10', (46) 

unstable particles and / ^ k then the probability theory states that the 
number Nsurv unstable particles 

Nsurv = V4t)\^^^^^Af^ ~ 10'-'^ > 0, (47) 
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has to survive up to time t ~ tas- Sources creating such numbers of unstable 
particles are known from cosmology and astrophysics. The Big Bang is the 
obvious example of such a source. Some other examples include processes 
taking place in galactic nuclei (galactic cores) and inside stars, etc. 

So let us assume that we have an astrophysical source creating a suffi- 
ciently large number of unstable particles in unit of time and emitting a flux 
of these particles and that this flux is constant or slowly varying in time. 
Consider as an example a flux of neutrons. From (!3T!) it follows that for the 
neutron f^^ ~ (250t„ — 300t„), where r„ ~ 886 [s]. If the energies of these 
neutrons are of order 30 x 10^^ [eV] then during time t ~ f^^ they can reach 
a distance rf" ~ 25000 [ly], that is the distance of about a half of the Milky 
Way radius. Now if in a unit ot time a suitably large number of neutrons jV„ 
of the energies mentioned is created by this source then in the distance dP' 
from the source a number of spherically symmetric space areas (halos) sur- 
rounding the source, where neutron instantaneous energies Sn{t) are much 
larger than its energy S^l = —^^^==, [m^ is the neutron rest mass and Vn 

denotes its velocity) have to appear (see Fig. ([8])). Of course this conclusion 
holds also for other unstable particles (/)„ produced by this source. 

Every kind of particles (pa has its own halos located at distances rf^° , 

dt-v^jt\ (fc=l,2,...), 

from the source. Radiuses d'^" of these halos are determined by the particles' 
velocities v^p^ and by times tf;,^ when instantaneous energies S^^{t) have local 
maxima. 

Unstable particles (pa forming these halos and having instantaneous en- 
ergies S^^{t) ^ £^9 = m9^ c^ have to interact gravitationally with objects 



outside of these halos as particles of masses m<^^ (t) = -^ £^^ (t) ^ m°^ . 
The possible observable effects depend on the astrophysical source of these 
particles considered. 

If the halos are formed by unstable particles emitted as a result of internal 
star processes then in the case of very young stars cosmic dust and gases 
should be attracted by these halos as a result of a gravity attraction. So, 
the halos should be a places where the dust and gases condensate. Thus in 
the case of very young stars one may consider the halos as the places where 
planets are born. On the other hand in the case of much older stars a presence 
of halos should manifest itself in tiny changes of velocities and accelerations 
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Figure 8: Halos surrounding a source of unstable particles 

of object moving in the considered planetary star system relating to those 
calculated without taking into account of the halos presence. 

If the halos are formed by unstable particles emitted by a galaxy core 

F (f) 

and these particles are such that the ratio -%^ is suitably large inside the 

halos, then rotational velocities of stars rounding the galaxy center outside 
the halos should differ from those calculated without taking into account the 
halos. Thus the halos may affect the form of rotation curves of galaxies. (Of 
course, we do do not assume that the sole factor affecting the form of the 
rotation curves are these halos). Another possible effect is that the velocities 
of particles crossing these galactic halos should slightly vary in time due to 
gravitational interactions, i. e. they should gain some acceleration. This 
should cause charged particles to emit electromagnetic radiation when they 
cross the halo. 

Note that the above mentioned effects seems to be possible to examine. 
All these effect are the simple consequence of the fact that the instantaneous 
energy S^f,^ (t) of unstable particles becomes large compared with £^^ and for 
some times even extremely large. On the other hand this property of S^^{t) 
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Figure 9: Halos generated by a few not too distant sources S^, S^. 

results from the rigorous analysis of properties of the quantum mechanical 
survival probability a{t) (see ([2]) ) and from the assumption that the energy 
spectrum is bounded from below. 
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